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\b" : ABSTRACT 

o 
o 

We investigate the dependences of galaxy star-formation history and galaxy 
morphology on environment, using color and Ha equivalent width as star- 
formation history indicators, using concentration and central surface brightness 
as morphological indicators, and using clustocentric distance as an environment 
indicator. Clustocentric distance has the virtue that it can be measured with 
very high precision over a large dynamic range. We find the following asymme- 
try between morphological and star-formation history parameters: star-formation 
history parameters relate directly to the clustocentric distance while morpholog- 
ical parameters relate to the clustocentric distance only indirectly through their 
relationships with star-formation history. This asymmetry has important impli- 
cations for the role that environment plays in shaping galaxy properties and it 
^ ■ places strong constraints on theoretical models of galaxy formation. Current 
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semi-analytic models do not reproduce this effect. 



1. Introduction 



It has been known for quite some time that the statistical properties of galaxies are 
closely related to their surrounding environments. Much of previous enviro nment-related 



work focused on the relationship between morpholo gy and environment (e.g., |Hubble 
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bulge-dominated galaxies are more strongly clustered than disk-dominated galaxies. Since 
spectroscopic and photometric properties of galaxies are strongly correlated with morphology, 



also functions of environment (e.g., 
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Since regions of the Universe with different densities evolve at different rates, we expect 
these environment dependencies to contain crucial information about galaxy formation and 
evolution. Considering that all of the observed statistical properties of galaxies appear to 
be interdependent, it is important to ask which properties are correlated with environment 
independently of the others; i.e., which environment relationships are fundamental and which 
ones are merely the products of other relations. 



Blanton et al.l (j2005al ) found that, out of color, luminosity, surface brightness, and radial 
concentration, the color and luminosity of a galaxy appear to be the only properties that are 
directly related to the local overdensity. Surface brightness and concentration appear to be 
related to the environment only through their relationships with color and luminosity. The 
main limitations in this study were caused by their environment indicator, local overdensity, 
since it has relative ly low signal-to-noise a n d it c annot resolve scales smaller than ~ 1 Mpc. 
In a similar study, IChristlein &: Zabludoffl (120051 ) used the clustocentric distance as an en- 
vironment indicator found that there is a residual dependence of ongoing star formation on 
environment among galaxies of similar morphology, stellar mass, and mean stellar age. This 
study was limited by a relatively small sample size (1,637 galaxies in 6 clusters) as well as 
by the star formation indicator, Oil, which is sensitive to dust and metallicity. 

In this shor t paper, we bu i ld upo n previous studies in several ways. We conduct a 
study similar to iBlanton et al.l (j2005al ). except that we use the clustocentric distance as 
our environment indicator. The clustocentric distance — the distance to the nearest galaxy 
cluster center — is a fundamental and precisely measurable environment indicator; indeed it 
can be measured at much highe r signal-to- nois e than e stimates of overdensity on fixed scales. 
Thus, we probe the findings of |R 
higher signal-to-noise. We use the 



anton et al. (2005a 



on much smaller scales and at much 
Berlind et al.l (120061 ) cluster catalog that was created with 



an algorithm that made no use of galaxy colors, so our galaxy colors are unbiased even in the 
cluster cores. The cluster catalog was designed to recover groups of galaxies that are bound 
in the same underlying dark matter halo. In order to determine whether our results can be 
easily understood within our current understanding of galaxy formation and evolution, we 
compare our observational results from the Sloan Digital Sky Sur yey to analogous re sults 
from the semi-analytic Millennium simulation of galaxy formation (ICroton et al.ll2006l ). 
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In what follows, a cosmological world model with (^Mj^a) 
Hubble constant is parameterized as H = 100 h kms _1 Mpc~ 
calculating distances and volumes h—1 except where otherwise noted (e.g.. lHogglll999f ) 



(0.3, 0.7) is adopted, the 
and for t he purpose s of 



2. Data and Analysis 



2.1. Observations 



The SDSS has taken uqriz CCD imaging and spectroscopy of many 1 5 galaxies at 
r < 17.77 mag fe.g- lGunn et alill998l: lYork et al.ll200ol : lstoughton et alJliooi l. All the data 



processing, i ncluding astrometr y (IPier et al.l 120031). source ident i fication, deblendin g and 



photometry OLupton et al 



troscopic target selection (lEisenstein et al.ll200ll; 



20011). calibration (|Fukugita et al 



1996 



Strauss et al. 



2002 



Smith et al.ll2002l), spec - 



Richards et al.M2002h 



spectroscopic fiber placement (IBlanton et al.ll2003al ). spectral data reduction and analysis 
(Schlegel & Buries, in preparation, Schlegel in preparation) are performed with automated 
SDSS software. 

Galaxy colors ar e computed in fixed bandpasses, using Galactic extincti on corrections 



(Sch 



egel et al.l Il998l ) and K corrections (computed with kcorrect v3_2; IBlanton et al. 



20031 ) . They are K corrected to the redshift z = observed bandpasses so they c an be directly 



comp ared to the rest-frame SDSS ugriz outputs of the Millennium simulation (ICroton et al. 
2006h . 



For the purposes of computing large-scale structure statistics, we have assembled a 
complete subsample of SPSS galaxies known as the NYU LSS samplel4. This subsample 
is described elsewhere ( IBlanton et al.l l2005bl ) ; it is selected to have a well-defined window 
function and magnitude limit. In addition, the galaxies in the sample used here were selected 
to have apparent magnitudes in the range 14.5 < r < 17.77 mag, redshift in the range 
0.015 < z < 0.068, and absolute magnitude Mo.ij > —24 mag. These cuts left 52,485 
galaxies. 

A seeing-convolved Sersic model is fit to the azimuthally aver aged radial profile of 
every galaxy in the observed-frame i band, as described elsewhere ( IBlanton et al.l l2003bl ; 
Strateva et al.l 1200 ll ). The Sersic model has surface brightness I related to angular radius r 
by / oc exp[— (r/ro)*- 1 /™- 1 ], so the parameter n (Sersic index) is a measure of radial concen- 
tration (seeing-corrected). At n = 1 the profile is exponential, and at n = 4 the profile is 
de Vaucouleurs. In the fits shown here, values in the range 0.5 < n < 5.5 were allowed. 



To every best-fit Sersic profile, the 



Petrosian 



(119761 ) photometry technique is applied, 
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with the same parameters as used in the SDSS survey. This supplies seeing-corrected Pet- 
rosian magnitudes and radii. A i^-corrected surface-brightness /io.o, in the ao z band is com- 
puted by dividing half the iT-corrected Petrosian light by the area of the Petrosian half-light 
circle. 

The Ha line flux is measured in a 20 A width interval centered on the line. Before 
the flux is computed , a best-fit model consisting of scaled old-galaxy and A-star spectra 



(jQuintero et al.l 120041 ) is scaled to have the same flux continuum as the data in the vicinity 
of the emission line and subtracted to leave a continuum-subtracted line spectrum. This 
method fairly accurately models the Ha absorption trough in the continuum, although in 
detail it leaves small negative residuals. The flux is converted to a rest-frame EW with 
a continuum found by taking the inverse-variance-weighted average of two sections of the 
spectrum about 150 A in size and on eit her side of the emission line. Further details are 



presented elsewhere (jQuintero et al.l 120041 ) . 



A caveat to this analysis is that the 3 arcsec diameter spectroscopic fibers of the SDSS 
spectrographs do not obtain all of each galaxy's light because at redshifts of 0.015 < z < 
0.068 they represent apertures of between 0.6 and 2.7 h^ 1 kpc diameter. Consequently the 
integrity of our Ha EW measurement is a function of galaxy size, inclination, and morphol- 
ogy. We have looked at variations in our Ha EW results as a function of redshift and found 
that the quantitative results differ but our qualitative results remain the same. 



We use the group and cluster catalog described in iBerlind et al.l (120061 ) . The catalog is 
obtained from a volume-limited sample of galaxies that is complete down to an r band 
absolute magnitude of M r < —19 mag and goes out to a redshift of 0.068. Groups are 
identified using a friends-of-friends algorithm (see e.g., Geller & Huchra 1983; Davis et al 
1985) with perpendicular and line-of-sight linking lengths equal to 0.14 and 0.75 times the 
mean inter-galaxy separation, respectively. These parameters were chosen with the help of 
mock galaxy catalogs to produce galaxy groups that most closely resemble galaxy systems 
that occupy the same dark matter halos. The resulting catalog contains 985 systems of 
richness iV ga i > 5 galaxies with an r-band absolute magnitude limit of M r < — 19 mag. For 
consistency, we call these objects "clusters". Note that, unlike some other catalogs, galaxy 
colors are not used in cluster identification. 

Berlind et al. (2006b, in preparation) calculate rough mass estimates for the clusters 
using the cluster luminosity function (where luminosity is defined as the total luminosity 
in M r < —19 mag galaxies in the cluster) and assuming a monotonic relation between a 
cluster's luminosity and the mass of its underlying dark matter halo. By matching the 
measured space density of clusters to the theoretical space density of dark matter halos 
(given the concordance cosmological model and a standard halo mass function), they assign 
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a virial halo mass to each cluster luminosity. The masses derived in this way ignore the 
scatter in mass at fixed cluster luminosity and are only meant to be approximate. The 
median mass estimate for the clusters used here is 3 x 10 13 solar masses. Each cluster has 
an associated "virial radius" of 



M 



^4tt200 p 0i 

where M is the estimated mass of the cluster and p Q is the current mean density of the 
Universe. By this method, the median estimated virial radius for the clusters used here is 
~ 1 Mpc. Note that this method for determination of the virial radii is very different from 
that employed by o ther investigators. Some have used a quasi-empi rical formula based on 



velocity dispersion (IGomez et al.ll2003l ; IChristlein &: Zabludofil 12005), oth ers have assumed 



that cluster mass is directly proportional to richness (ILewis et al. 



20021 ) ; in general these 



methods differ substantially, and produce cluster catalogs with very different mass functions. 



We use the cluster centers given by iBerlind et al.l (120061 ) , which are computed as the 
mean of the member galaxy positions. We then calculate the transverse projected clustocen- 
tric distance D c \ from each galaxy to its nearest cluster center on the sky within ±1000 km s _1 
in radial velocity. 



2.2. Simulations 

To better understand our observations, we do a similar analysis for the Millennium sim- 
ulation galaxy catalog. The Millennium simulation is a large cosmological N-body simulation 
of 2160 3 dark matter particles in a 500 3 /i~ 3 Mpc 3 box. The formation of galaxies is simu- 
lated by using a semi-analytical model whose parameters and assumptions are tweaked in 
order to reproduce the joint luminosity-color, morphology, gas mass, and cent ral black hole 



mass distributions of low-redshift galaxies. For a more detailed description, see lCroton et al. 



fl2006f ) 



The publicly available output of this simulation (ICroton et al.ll2006l ) is a catalog of ~ 9 
million galaxies that contains positions, velocities, bulge and total SDSS-band magnitudes, 
bulge and total masses, gas masses (cold, hot, and ejected), black hole masses, and star 
formation rates. 

From this catalog we use the ratio of bulge-to-total luminosity as a morphology indi- 
cator. This B/T value is an analog to our observed concentration and surface brightness. 
Similarly, we use the star formation rate SFR and SDSS-band °'°[g — r] color from this cata- 
log as star- formation history indicators; these are analogs to our Ha EW and 0,0 [g — r] color 
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measurements, respectively. 



Before running the cluster- finding algorithm (jBerlind et al.ll2006l ). we distort the Mil- 
lennium coordinates to mimic redshift-space distortions by changing the output coordinates 
of each galaxy from [(x), (y), (z)] to [(x), (y), (z + jf^)}- We then run the same friends-of- 
friends algorithm that we did on the data, using the perpendicular linking length in the x— 
and y— directions and the line-of-sight linking length in the z-direction. After creating this 
cluster catalog we calculate each galaxy's clustocentric distance using a similar algorithm we 
use on the data: We calculate the transverse projected clustocentric distance D d from each 
galaxy to its nearest cluster center on the x — y plane within ±10Mpc in the z-direction. 



3. Results 

Figure [T] shows how radial concentration and color relate to clustocentric distance within 
narrow subsamples of each other. For the 2x6 grid on the left, the top- left panel shows how 
the 5, 25, 50 (in bold), 75, and 95 percentiles of the concentration distribution vary with 
clustocentric distance for the entire sample. The color distribution for the entire sample is 
shown in the right panel. The gradient in the concentration percentiles show that galaxies 
with small clustocentric distances tend to be more concentrated. The subsequent rows show 
the same relation but for narrow color subsamples whose color distribution is shown in the 
right panels. Notice the gradient seen in the top panel nearly vanishes for the narrow color 
subsamples below. The absence of a concentration gradient for these subsamples shows 
that at fixed color, there is no dependence of concentration on environment. The right 
2x6 grid shows the same plot but with concentration and color interchanged. Notice that 
a color gradient remains when looking at narrow concentration subsamples. This shows 
that at fixed concentration, there is a residual dependence of color on environment. The 
asymmetry seen in this figure is also seen in Figure [2] which uses surface brightness, instead 
of concentration, as a morphology tracer. This as ymmetry is the princip al result of this 



paper. This asymmetry is similar to that found by iBlanton et al.l (l2005al ). but it appears 



here at much higher significance and over a larger dynamic range. 

Figure [3] also shows the same asymmetry shown in Figure [1] but using Ha EW, instead 
of color, as our star-formation-history param eter. Note that the right 2x6 grid agrees with 



the findings of IChristlein fc Zabludoffl (120051 1 . 



Figure H] displays the asymmetry shown in Figure [T] in a different way. The left panel 
shows that out of color and clustocentric distance, concentration depends only on color since 
the mean concentration contours are perpendicular to the color axis. The center panel shows 
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that color depends on both concentration and clustocentric distance. The right panel shows 
that clustocentric distance depends more on color than concentration, since the contours are 
nearly perpendicular to the color axis. 

Figure [5] is analogous to Figure [I] but showing the redshift zero outputs from the Mil- 
lennium simulation. Instead of concentration n we show the ratio of bulge luminosity to 
the total luminosity B/T. Notice that the asymmetry seen in the observed data is not 
reproduced here. The asymetry is also absent when we replace color with star-formation 
rate. 

Figure [6] is the same as Figure [5] but with 30 percent random errors added to the B/T 
values. The "random errors" were produced by randomly adding a Gaussian distribution 
with a dispursion of a = 0.3 to the Millennium B/T values. Notice that the behavior of 
the trends do not significantly change when adding these large errors to B/T. This figure 
suggests that the asymmetry observed in the SDSS data is not due to the fact that we 
measure star-formation histories with much higher precision than we measure morphologies. 

Another interesting result from Figures 1-3 is that the dependence of star-formation- 
history and morphological parameters on clustocentric distance is mostly limited to within 
the cluster virial radius. The top-left panels of these figures show that the trends observed at 
distances of less than the virial radius (denoted by the vertical dashed lines) nearly disappear 
at larger radii, where the relations between star- formation/morphological parameters and 
clustocentric distance are mostly flat. This result is not reproduced by the Millenium semi- 
analytic model. The top-left panels of Figure [5] show that, though the relations between 
model galaxy parameters and clustocentric radius flatten out near the virial radius, the 
strong trends resume at larger radii. 



4. Summary and Discussion 

Using a complete sample of 52,485 galaxies in the redshift range of 0.015 < z < 0.068 we 
examine how star- formation history parameters — [g — r] color and Ha EW — and morphol- 
ogy parameters — radial concentration n and surface brightness /io.Oj — relate to environment. 
We use the transverse projected distance to the nearest cluster— the clu stocentric distance 



D c i — as our environment parameter. Similar to iBlanton et al.l (j2005al ). we find that our 
morphology tracers (concentration and surface brightness) appear to be related to environ- 
ment (D c \ in this case) only indirectly through their relationships with the star-formation 
history tracers. This suggests that the well known morphology-environment relation is a 
result of the star- formation-history-environment relation. We note that, although our sim- 
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ple morphology indicators do not show direct dependence on clustocentric distance, this 



result does not disagree wi t h that the original mor phology-density studies (jHubbld Il936 



Oemlerlll974j ; lDresslerlll980l ; iPostman fc Gellerlll984l ) because these studies did not attempt 
to separate morphol ogical and star-fo rmation dependences. Our results also do not disagree 
with previous work (IVogt et al.l I2004J ) which suggests that gas stripping plays a significant 
role in the morphological transformation and rapid truncation of star formation by show- 
ing asymmetries in HI and Ha flux on the leading edge of infalling spiral galaxies. There 
is no disagreement because our very blunt morphology indicators are insensitive to those 
asymmetries, and the asymmetries are only visible in a small fraction of all cluster galaxies. 



Where our resul ts overlap those of previous investigators (IChristlein &: Zabludofj 12005 



Blanton et al.ll2005al ). we mostly find good agreement. In particular, we find that galaxy 
properties depend on clu stocentric distance in much the same way as they do on other 
environmental indicators (IBlanton et al.ll2005al ). Our result improves on this previous one 
because the clustocentric distance is measured at much higher signal-to-noise and probes 
environments on scales m uch smaller than local o yerde nsity measurements. Our results 
improve on the findings of IChristlein fc Zabludofll ((20051) by reproducing their result with 
far higher signal to noise and by also by examining the inverse of their study. They found 
that star formation rate is related to environment independently of morphology. We find 
that both color and star formation rate behave in this way as well as that morphology doesn't 
have a dependence on environment independent of star-formation history. 

It is worth noting that this asymmetry may, in principle, be attributed to the mea- 
surements that were made instead of the properties of galaxies. This is possible because 
our star-formation-history properties are measured with far higher precision than our mor- 
phology properties. However, we do not believe this is true for two main reasons; (1) when 
we add large asymmetric noise to the Millennium bulge-to-total luminosities we still find 
no clear asymmetry and (2) the fact that we observe this asymmetry regardless of which 
star-formation history and morphology tracers we use, strongly suggests that it is a descrip- 
tion of galaxy properties themselves. We stress that our main result, that there is no direct 
morphology-environment relation, holds for morphology defined by concentration and sur- 
face brightness. It remains untested whether other characterizations of morphology, such as 
spiral structure or asymmetry in the shape of the light distribution, also show this trait. 

Our results also show that, where there are trends between galaxy properties and clus- 
tocentric distance, these trends are primarily restricted to within the cluster virial radii. 
This suggests that environmental processes that affect galaxy properties are fairly local (i.e., 
within the dark matter halo) and do not "act at a distance ". These results are not in 



agreement with those of iGomez et al.l (120031 ) and iLewis et al.l (120021 ) who found that there 
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are trends out to ~ 3 times the virial radius. The discrepancy can be attributed to our 
different ways of estimating cluster virial radii. These authors use the iGirardi et al.l (119981 ) 
appro ximation relating viri al radius to velocity dispersion, whereas we use virial radii based 
on the Berlind et al. I d2006h ma s s estim ates, which are derived from cluster abundances. We 
have checked the IGirardi et al.l (119981 ) approximation on our clusters and find that it yields 
virial radii that are ~ 2 times smaller than the ones we use, which explains much of the 
discrepancy. We trust our mass estimates bec ause they are based on abundances, which 
have been shown to be unbiased for our clusters iBerlind et al.l (120061 ). Masses based on the 
Girardi et al.l (119981 ) approximation yield an incorrect cluster mass function. 



The asymmetry that we see between star-formation and morphological parameters has 
interesting implications for the role that environment plays in shaping galaxy properties. 
Our results suggest that environmental processes can affect galaxy star-formation histories 
without simultaneously affecting morphology, but in cases where an environmental mech- 
anism does affect galaxy morphology, it must also affect star-formation history. Another 
possible explanation is that environment affects both star-formation history and morphol- 
ogy, but the changes in star-formation history occur on a much more rapid timescale than 
the corresponding changes to morphology. 

In this short paper, we have investigated how galaxy star-formation histories and mor- 
phologies depend on environment. We find the following asymmetry: At fixed color, there is 
no dependence of morphology on environment, while at fixed morphology, there is a residual 
dependence of color on environment. This result suggests that the morphology-environment 
relation is a result of the combination of the star-formation-history-environment and the 
morphology-star-formation-history relations. Moreover, we performed a similar analysis on 
a semi-analytic model of galaxy formation and find that this asymmetry is not reproduced. 
Our results have important implications for the role that environment plays in shaping 
galaxy star-formation histories and morphology, and they place strong constraints on mod- 
els of galaxy formation. 
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Fig. 1. — The dependence of radial concentration n (left grid) and 00 [g — r] color (right grid) 
on clustocentric distance D c \. For the 2x6 grid on the left, the top-left panel shows how the 
5, 25, 50 (in bold), 75, and 95 percentiles of the distribution of concentration n (Sersic index) 
depend on clustocentric distance D c \ for the entire sample. The top-right panel of this grid 
shows the °-°[g — r] color distribution of the entire sample. The errors in the 50 percentile, 
calculated using the jackknife method (10 jackknife trials in each of which a contiguous 1/10 
of the survey footprint is dropped) are overplotted in white (they are tiny). The subsequent 
rows show the same but for narrow color subsamples (color quintiles). The 2x6 grid on 
the right is very similar but with the n and 00 [g — r] properties interchanged. The left 
2x6 grid shows that in color subsamples, there is no residual dependence of concentration 
on environment, while the right 2x6 grid shows that in concentration subsamples, there is 
a residual dependence of color on environment. This asymmetry is the principal result of 
this paper. 
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Fig. 2. — Similar to Figure HJ except replacing the radial concentration n with surface 
brightness /zo.Oj. Notice the asymmetry seen in Figure [1] also occurs here. 




Fig. 3. — Similar to Figured! except replacing the °'° [g—r] color with our other star- formation 
history indicator: Ha emission line equivalent width (Ha EW). Notice the asymmetry seen 
in Figure [1] also occurs here. 
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Fig. 4. — The left panel shows contours of mean radial concentration n as a function of 
°-°[g — r] color and clustocentric distance D c \. The mean is computed in a sliding box shown 
in the center of the plot. The center panel shows contours of mean color as a function of 
concentration and clustocentric distance. The right panel shows contours of mean clusto- 
centric distance as a function of color and concentration. Note that the contours in the left 
and right panels are more horizontal than the contours in the center panel. 



-17- 




Fig. 5. — The same 2x6 grids shown in Figured] but showing the Millennium simulation 
redshift zero outputs. Instead of concentration n we show the ratio B/T of the bulge lumi- 
nosity to total luminosity as a morphology indicator. Notice the asymmetry seen in Figured] 
is not seen here. 
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Fig. 6. — The same as Figure [5] but with thirty percent scatter (see text) added to the 
bulge-to-total values. Notice that the behavior of the trends do not significantly change 
when adding these large errors to B/T. 



